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A new parametrization PC-PK1 for the nuclear covariant energy density functional with nonlin- 
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illustrated in the description for infinite nuclear matter and finite nuclei including the ground-state 
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dence of binding energy along either the isotopic or the isotonic chains, which makes it reliable 
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are determined by constrained calculations for triaxial shapes. 
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I. INTRODUCTION 



In the past years, the unstable nuclear beams have extended our knowledge of nuclear 
physics from the stable nuclei to the unstable nuclei far from the stability line — so-called 
"exotic nuclei" . Extensive research in this area shows a lot of entirely unexpected features 
and novel aspects of nuclear structure such as the halo phenomenon [1-3], the disappearance 
of traditional magic numbers and the occurrence of new ones [4]. The exotic nuclei play 
important roles in nuclear astrophysics, since their properties are crucial to stellar nucle- 
osynthesis. To understand the physics in exotic nuclei, it becomes very important to find a 
reliable theory and improve the reliability for predicting the properties of more exotic nuclei 
close to proton and neutron drip lines. 

Nuclear energy density functional (EDF) theory [5] has played an important role in a 
self-consistent description of nuclei. With a few parameters, EDF theory is able to give a 
satisfactory description for the ground state properties of spherical and deformed nuclei all 
over the nuclide chart. Detailed discussion on the EDF theory can be seen in Ref. [6] for 
nonrelativistic representations and in Refs. [7, 8] for relativistic ones. 

There exist a number of attractive features in the covariant EDF theory, especially in 
its practical applications of self-consistent relativistic mean-field (RMF) framework [7, 8]. 
The most obvious one is the natural inclusion of the nucleon spin degree of freedom and the 
resulting nuclear spin-orbit potential that emerges automatically with the empirical strength 
in a covariant way. The relativistic effects are responsible for the empirical existence of 
approximate pseudospin symmetry in the nuclear single-particle spectra [9]. Moreover, a 
covariant treatment of nuclear matter provides a distinction between scalar and four-vector 
nucleon self energies, leading to a natural saturation mechanism. 

The most widely used RMF framework is based on the finite-range meson-exchange repre- 
sentation (RMF-FR), in which the nucleus is described as a system of Dirac nucleons which 
interact with each other via the exchange of mesons. The isoscalar-scalar a meson, the 
isoscalar-vector u meson, and the isovector-vector p meson build the minimal set of meson 
fields that, together with the electromagnetic field, is necessary for a description of bulk 
and single-particle nuclear properties. Moreover, a quantitative treatment of nuclear mat- 
ter and finite nuclei needs a medium dependence of effective mean-field interactions, which 
can be introduced by including nonlinear meson self-interaction terms in the Lagrangian 



2 



or by assuming explicit density dependence for the meson-nucleon couplings. Of course, at 
the energy characteristic for nuclear binding and low-lying excited states, the heavy-meson 
exchange (cr, u, p) is just a convenient representation of the effective nuclear interaction. 

Since the exchange of heavy mesons is associated with short-distance dynamics that can- 
not be resolved at low energies, as an alternative, the relativistic point-coupling (RMF-PC) 
model [10, 11] is proposed by using the zero-range point-coupling interaction instead of the 
meson exchange, i.e., in each channel (scalar- isoscalar, vector- isoscalar, scalar- isovector, and 
vector-isovector) meson exchange is replaced by the corresponding local four-point (contact) 
interaction between nucleons. Analogously, in the case of contact interactions, the medium 
effects can be taken into account by including higher-order (nonlinear coupling) interac- 
tion terms or by assuming a density dependence of strength parameters for the coupling 
interactions. 

In recent years, the RMF-PC model has attracted more and more attentions due to 
the following advantages. Firstly, it avoids the possible physical constrains introduced by 
explicit usage of the Klein-Gordon equation to describe mean meson fields, especially the 
fictitious a meson. Secondly, it is possible to study the role of naturalness [12, 13] in effective 
theories for nuclear structure related problems. Thirdly, it provides more opportunities to 
investigate its relationship to the nonrelativistic approaches [14] . Finally, it is relatively easy 
to study the effects beyond mean-field for the nuclear low-lying collective excited states. 

In practical application of the RMF-PC model, the most widely used nonlinear cou- 
pling parameterizations include PC-LA [10] and PC-F1 [11]. PC-LA is determined by the 
ground-state observables of 16 0, 88 Sr, and 208 Pb. Due to the explicit omission of the pairing 
interaction, the pairing effects are not included in the fitting procedure. Moreover, the test 
for naturalness in Ref. [12] shows that only six of the nine coupling constants are natural. As 
an improvement, PC-F1 is optimized to observables of 17 spherical nuclei including open- 
shell nuclei, and the pairing correlation is considered through a standard BCS approach in 
the fitting procedure. Furthermore, all the coupling constants of PC-F1 are turned out to 
be natural [11]. However, the isospin dependence of binding energy given by PC-F1 along 
either the isotopic or the isotonic chains deviates from the data remarkably. 

Recently, a density-dependent parametrization DD-PC1 is proposed from the equation 
of state (EOS) of nuclear matter and the masses of 64 axially deformed nuclei in the mass 
regions A ~ 150 — 180 and A ~ 230 — 250 [15]. Although it reproduces the binding energies, 
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deformations, and charge radii of deformed nuclei quite well, the differences between the 
predicted binding energies and the corresponding data are somewhat large for spherical 
nuclei. 

Therefore, it is necessary to have a new parametrization for the nuclear covariant energy 
density functional with point-coupling interaction to describe both the nuclear matter and 
finite nuclei properties. In this work, a new parametrization PC-PK1 with nonlinear coupling 
interactions is proposed. In Sec. II, the theoretical framework for the relativistic point- 
coupling model is briefly outlined. The numerical details are given in Sec. III. In Sec. 
IV- VII, a series of illustrative descriptions for the nuclear matter, spherical nuclei, deformed 
nuclei as well as the nuclear excited properties are presented. Finally, a summary is given 
in Sec. VIII. 

II. THEORETICAL FRAMEWORK 

The basic building blocks of RMF theory with point-coupling vertices are 

(ipOTip), Oe{l,r}, Te {1,7^,75,757^,^}, (1) 

where ip is Dirac spinor field of nucleon, f is the isospin Pauli matrix, and V generally 
denotes the 4x4 Dirac matrices. There are ten such building blocks characterized by their 
transformation characteristics in isospin and Minkowski space. In this paper, vectors in the 
isospin space are denoted by arrows and the space vectors by bold type. Greek indices /i 
and v run over the Minkowski indices 0, 1, 2, and 3. 

A general effective Lagrangian can be written as a power series in ifjOTip and their 
derivatives. We start with the following Lagrangian density of the point-coupling model 

£ = £ frcc + £ 4f + £ hot + £ dcr + £ cm , (2) 

which is divided as the Lagrangian density for free nucleons £ free , 

£ free = i>(ij^ - m)1>, (3) 

the four-fermion point-coupling terms £ 4f , 

1 - - 1 - 

£ = -^sWjfW) - 2 a ^(V'7W')(V'7 / V) 

--«TS'(^T'0)(^f^) - 2 a rv , (^7 M V0(^7 M V')j ( 4 ) 
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the higher order terms £ which are responsible for the effects of medium dependence, 

£ hot = -iflKW - \is{^Y - \iv[$iMh^)\\ (5) 

the gradient terms £ der which are included to simulate the effects of finite-range, 

~8 TS d v {^)d v {^) - ^ TV ^(^f 7 ^)^(^7^), (6) 

and the electromagnetic interaction terms £ em , 

£ cm = -\f^F, u - et-^fjAr (7) 

For the Lagrangian density in Eq. (2), m is the nucleon mass and e is the charge unit for 
protons. and F^ u are respectively the four-vector potential and field strength tensor of 
the electromagnetic field. There are totally 11 coupling constants, as, ay, «ts, Q-tv, Ps, Is, 
7v, 5s, 5y, 5ts, an d 5tv, m which a refers to the four-fermion term, (3 and 7 respectively the 
third- and fourth-order terms, and 5 the derivative couplings. The subscripts S, V, and T 
respectively indicate the symmetries of the couplings, i.e., S stands for scalar, V for vector, 
and T for isovector. 

From former experience [11], we neglect the isovector-scalar channel in Eq. (2) since a 
fit including the isovector-scalar interaction does not improve the description of nuclear 
ground-state properties. Consequently, there are nine free parameters in the present RMF- 
PC model, which are comparable with those in the RMF-FR model. Furthermore, the 
pseudoscalar 75 and pseudovector 757^ channels are also neglected in Eq. (2) since they do 
not contribute at the Hartree level due to parity conservation in nuclei. 

Similar to the RMF-FR case, the mean-field approximation leads to the replacement of 
the operators ^(OT)^ in Eq. (2) by their expectation values which become bilinear forms 
of the nucleon Dirac spinor ip^, 

^{OY)^ (<J#(OT)^|$) = vlMOThiJk, (8) 

k 

where i indicates S, V, and TV. The sum runs over only positive energy states with the 
occupation probabilities v\, i.e., the "no-sea" approximation. Based on these approxima- 
tions, one finds the energy density functional for a nuclear system 

*W[t,/>5X4J = / d 3 rS(r), (9) 
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with the energy density 

S(r) = £ kh \r) + £ int (r) + £ cm (r), (10) 
which is composed of a kinetic part 

£kin (r) = J- v\ 4(r) (a • p + /3m) ^(r), (11) 

an interaction part 

+fw" + f(w") 2 + |^A J " (12) 

+—Jtv ■ Utv)h + — 3tv ■ A Utv)^ 
with the local densities and currents 

Ps(r) = ^2v 2 k tp k (r)^ k (r), (13a) 

k 

tt(r) = J2 v lMr)YMr), (13b) 

k 

JTv(r) = J2 v lMr)r^Mr), (13c) 

k 

and an electromagnetic part 



£ cm (r) = -F^ - F°»d A, + eA^. (14) 

Minimizing the energy density functional Eq. (9) with respect to ipk, one obtains the 
Dirac equation for the single nucleons 

[ lli (id»-V i )-(m + S)]if> k = 0. (15) 

The single-particle effective Hamiltonian contains local scalar S(r) and vector V M (r) poten- 
tials, 

S(r) = T s , V^r) = T^ + r-T%, (16) 

where the nucleon scalar-isoscalar Tg, vector-isoscalar T^ 1 , and vector-isovector T^ v self- 
energies are given in terms of the various densities, 

Ss = agpg + [3gpg + 'Jgpg + SgApg, (17a) 

T» = a V 3» + lv (Jv) 3 + 8v&& + eA", (17b) 
T^ v = a TV j£ v + 6 TV Aj^y. (17c) 



For a system with time reversal invariance, the space-like components of the currents ji in 
Eq. (13) and the vector potential V(r) in Eq. (16) vanish. Furthermore, one can assume 
that the nucleon single-particle states do not mix isospin, i.e., the single-particle states are 
eigenstates of r 3 . Therefore only the third component of isovector potentials ^ TV survives. 
The Coulomb field A Q is determined by Poisson's equation. 

In addition to the self-consistent mean-field potentials, for open-shell nuclei, pairing cor- 
relations are taken into account by the BCS method with a smooth cutoff factor f k to 
simulate the effects of finite-range [16, 17], i.e., we have to add to the functional Eq. (9) a 
pairing energy term of the form depending on the pairing tensor k, 

E pair [K,K*} = fkfk'(kk\V PP \k'k')KlK k ,, (18) 

kk'>0 

with the smooth cut-off weight factor 

fk = l + exp[(e k -e F -AE T )/f, T y (19) 
where e k is the eigenvalue of the self-consistent single-particle field, and tp is the chemical 
potential determined by the particle number, ($|iV T |$) = iV r , with N T the particle number 
of neutron or proton. The cut-off parameters AE T and \x T = AE T /10 are chosen in such a 
way that 2 ^ f k = N T + 1.65iV r 2/3 [17]. 

fc>0 

In the following calculations, a density-independent 5-force in the pairing channel is 
adopted. Thus, the pairing energy is given by 

E paiT [K,K*} = ~J2^f d 3 rK* T (r)K T (r), (20) 

r=n,p 

where V T is the constant pairing strength and the pairing tensor n(r) reads 

K{r) = -2Y,fkU k v k \Mr)\ 2 - (21) 

fc>0 

The pairing strength parameters V T can be adjusted by fitting the average single-particle 
pairing gap 

(A) = Efc fk u k"Vk^k (22) 
Efc fkU k V k 

to the data obtained with a five-point formula. 

As the translational symmetry is broken in the mean-field approximation, proper treat- 
ment of center-of-mass (cm.) motion is very important and here the cm. correction energy 



is calculated by microscopic cm. correction 

B * = -aSA<^>- < 23) 

with A mass number and P c . m . = Ylf Pi ^ ne total momentum in the cm. frame. It has been 
shown that the microscopic cm. correction provides more reasonable and reliable results 
than phenomenological ones [18-20]. 

Therefore, the total energy for the nuclear system becomes 

E tot = E DF [r,p s ,j^A^ + E pah [K,K*] + E^, (24) 



III. NUMERICAL DETAILS 



In this work, a series of calculations have been performed for both the spherical and 
deformed nuclei. The Dirac equation for nucleons is solved in a three-dimensional harmonic 
oscillator basis [21]. For spherical calculations, by increasing the fermionic shells from Nf = 
20 to Nf = 22, the binding energy, charge radius, and neutron skin thickness in 208 Pb 
change by 0.003%, 0.007%, and 0.1% respectively. For 240 U, the binding energy changes 
by 0.001% from the axially deformed calculation with Nf = 16 to Nf = 18. Therefore, a 
basis of 20 major oscillator shells is used in the spherical calculations and 16 shells in the 
axially deformed cases. The triaxial calculations are performed with Nf = 12, which, for 
Nd isotopes, provides an accuracy of 0.04% for binding energies in comparison with the 
calculations with Nf = 14. To achieve an accuracy of ~ 100 keV in the description of both 
the fission barrier and the energy of the isomer state in 240 Pu, a basis of 20 oscillator shells 
has been adopted in both the axial and triaxial calculations. 

In order to determine the parameters of Lagrangian density in Eq. (2) and the pairing 
strength in Eq. (20), a multiparameter fitting to both the binding energies and charge radii 
for selected spherical nuclei is performed with the Levenberg-Marquardt method [22]. As 
usual, the masses of neutron and proton are fixed as 939 MeV. The corresponding data [23- 
25] for selected spherical nuclei used in the fitting procedure are listed in Table II and III. 
The empirical neutron pairing gaps for 122 Sn, 124 Sn, and 200 Pb as well as the proton ones 
for 92 Mo, 136 Xe, and 144 Sm obtained with five-point formula are also employed to constrain 
the pairing strengths. 



8 



With the experimental observable 0^ xp and the calculated value O™ 1 , by minimizing the 
square deviation 



r C> exp - 0, cal fa) 



X 2 (a) 



(25) 

the ensemble of parameters a can be obtained. Furthermore, in order to balance the influence 
of different observables, the weight is introduced for binding energies, charge radii, and 
empirical pairing gaps respectively. The corresponding weight Ui is roughly determined by 
the desired accuracy. Here the weights Ui are respectively 1.00 MeV for binding energies, 
0.02 fm for charge radii, and 0.05 MeV for empirical pairing gaps. A new parameter set 
PC-PK1, which contains the nine coupling constants in Eq. (2) and the pairing strength in 
Eq. (20), is obtained and listed in Table I. 

By scaling the coupling constants in accordance with the QCD-based Lagrangian, the 
naturalness in effective theories can be investigated [12, 13]. According to the QCD-based 
Lagrangian [13], 



C ~ -Ci r 



-i I r 



X 



(26) 



with ip the nucleon field, f n = 92.5 MeV the pion decay constant, and A = 770 MeV a generic 
QCD large-mass scale respectively, by taking into account the role of chiral symmetry in 
weakening iV-body forces byA = / + ra — 2^0 [26, 27], it has been found that six of the nine 
coupling constants in PC-LA and all of them in PC-F1 are natural, i.e., the QCD-scaled 
coupling constants q n are of order unity [11, 12]. 

Similarly, the nine coupling constants of PC-PK1 are also tested for the naturalness and 
all the dimensionless coefficients Q n are of order 1, as shown in the last column of Table I, 
which indicates that all the coupling constants in PC-PK1 are natural. 

Tables II and III list respectively the binding energies and charge radii for nuclei selected 
in the determination of PC-PK1, PC-F1, PC-LA, and NL3* [28] effective interactions. The 
corresponding root mean square (rms) deviation A together with the root of relative square 
(rrs) deviation 5 for the binding energy and charge radius are given in the last two rows 
of Table II and III, respectively. Compared with the other effective interactions, the newly 
obtained PC-PK1 provides a much better description for the experimental binding energies 
and the same good description for the charge radii. 
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IV. NUCLEAR MATTER PROPERTIES 



In this section, we will present the saturation properties and the equation of state (EOS) 
for nuclear matter in the covariant EDF with PC-PK1. The results will be compared with 
the corresponding empirical values as well as the predictions with PC-LA, PC-F1, DD-PC1, 
NL3*, and PK1 [19]. 

A. Saturation properties 

The saturation properties, including the binding energy per nucleon E/A, saturation 
density po, incompressibility Kq, nucleon effective mass M|> and M£, symmetry energy 
Egym, as well as the characteristics L and -K" asy for the density dependence of E sym will be 
investigated. 

There are several kinds of nucleon effective mass [29, 30]. Here we mainly focus on the 
Dirac mass and Landau mass. The Dirac mass M* D is defined through the nucleon scalar 
self-energy in the Dirac equation, i.e., = M + S5. It is directly related to the spin-orbit 
potential in finite nuclei and is thus a genuine relativistic quantity without nonrelativistic 
correspondence. While the Landau mass M£ = pdp/dE is related to the density of state both 
in relativistic and non-relativistic models. In relativistic models, the relation between the 
Dirac mass and the Landau mass is M£ = \Jpp + Mff , where pp is the Fermi momentum. 

The density dependence of the nuclear symmetry energy is very important to understand 
the properties of exotic nuclei with extreme isospin values, in particular the slope L = 
3po(dE sym /dp) p=po and curvature K sym = 9p 2 l (d 2 E S y Ta /d 2 p) p=Po of the symmetry energy at 
the saturation density p . In Refs. [31, 32], the isospin-dependent part, A^y « K sym — 6L, 
in the isobaric incompressibility K(5) = K + A asy 5 2 (with 5 = (p n — p p )/p), is often used 
to characterize the density dependence of the symmetry energy as both L and A asy can be 
extracted from the experiment empirically (see Ref. [33] and references therein). 

In Table IV, the saturation properties for nuclear matter, including the binding energy 
per nucleon E/A, saturation density po, incompressibility Kq, nucleon effective mass M£> 
and M£, symmetry energy E sym , as well as the characteristics L and K^y for the density 
dependence of E sym , predicted by PC-PK1 are listed in comparison with those by both 
point-coupling DD-PC1, PC-F1, PC-LA and meson exchange NL3*, PK1 sets. In general, 



10 



PC-PK1 gives good description for the saturation properties of nuclear matter. In particular, 
the predicted values for binding energy per nucleon and density at the saturation point are 
—16.12 MeV and 0.154 fm~ 3 , which agree well with the empirical values — 16 ± 1 MeV and 
0.166 ± 0.018 fm -3 [34], respectively. Moreover, the incompressibility given by PC-PK1 is 
238 MeV. 

For the effective masses, all the effective interactions give reasonable values between 
0.55 and 0.60 for the Dirac mass M^/M [35] as required by the spin-orbit splitting data 
in finite nuclei, but smaller Landau masses M* L jM compared with the empirical constraint 
0.8 ±0.1 [36], implying that they would give a small single-particle level density at the Fermi 
energy in finite nuclei as compared with data. 

The symmetry energies in the calculations with the nonlinear effective interactions PC- 
Fl, PC-LA, NL3* and PK1 are always larger than the empirical value (around 32 MeV) by 
around 16-18%, which is reduced to 11% for PC-PK1. Therefore all the interactions would 
predict large neutron skin thicknesses in finite nuclei except DD-PC1, which is adjusted by 
fixing E sym = 33 MeV. Moreover, the empirical L (88 ± 25 MeV) [33] and K &SJ (-550 ± 
100 MeV) [37] have been reproduced quite well by PC-PK1. 

B. Equation of state 

In Fig. 1, the binding energy per nucleon E/A for nuclear matter as a function of the 
baryon density pB given by PC-PK1 is shown in comparison with those by DD-PC1, PC-F1, 
PC-LA, NL3*, and PK1. All the effective interactions predict the similar E/A behavior 
with density below ps = 0.20 fm -3 due to the constraints from the properties of the finite 
nuclei. Divergence appears at supra-saturation densities, especially for the results given by 
PC-LA. It implies that the properties of finite nuclei are not sufficient in the determination 
of EDF to describe the EOS at supra-saturation densities that are directly related to the 
maximal mass of neutron star. The prediction by PC-PK1 is consistent with those by 
PK1 and PC-F1, while softer than those by NL3*. The ab-initio variational calculation 
for the symmetric nuclear matter [38] is also given for comparison, which coincides with 
the relativistic EOS with density below ps = 0.20 fm~ 3 but predicts softer EOS behavior 
at supra-saturation densities than the relativistic ones, except those given by PC-LA and 
DD-PC1. One should note that DD-PC1 has been adjusted to the EOS given by ab-initio 
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variational calculations [15]. 



V. SPHERICAL NUCLEI 

In this section, we will present the binding energies, two-neutron separation energies, 
single-particle levels, charge radii and neutron skin thicknesses for selected spherical isotopes 
and isotones in different mass regions in the covariant EDF with PC-PK1. The results will 
be compared with the corresponding data available as well as the predictions with DD-PC1, 
PC-F1, PC-LA, and NL3* sets. 

A. Binding energy 

The binding energies for the Ca, Ni, Sn, and Pb isotopes are calculated with PC-PK1 
and their deviations from the data [23] are shown in Fig. 2 in comparison with those with 
DD-PC1, PC-F1, PC-LA, and NL3*. The calculation with PC-PK1 reproduces the experi- 
mental binding energies within 1 MeV for the Ca isotopes and 2 MeV for both Sn and Pb 
isotopes. For Ni isotopes, although remarkable improvement is achieved by PC-PK1 in com- 
parison with results by the other interactions, there is still an underestimation of 2.4 — 5.2 
MeV for the even-even 58_64 Ni. Former investigations have shown that these nuclei are soft 
against deformation [39]. Therefore, the dynamic correlation energies gained by restoration 
of rotational symmetry and configuration mixing are expected to reduce these deviations. 
Similarly, the underestimations of the binding energies for neutron-deficient Pb isotopes can 
also be improved by configuration mixing, as demonstrated in the non-relativistic calcula- 
tions for the even-even 182 - 194 pb [40]. For the spherical Ca and Sn isotopes, the energies 
gained from the restoration of rotational symmetry and configuration mixing are expected to 
be much smaller as illustrated in the systematic beyond mean-field studies [41, 42]. There- 
fore, the inclusion of these energies will not change significantly the discrepancy between 
the mean- field results and the corresponding data for such spherical isotopes.. 

The binding energies for the isotonic chains are very important to examine the balance 
between the Coulomb field and the isovector channel of the Lagrangian density in Eq. (2). 
Here the binding energies for the N = 20, N = 50, N = 82, and N = 126 isotones are 
calculated with PC-PK1 and their deviations from the data [23] are shown in Fig. 3 in com- 
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parison with those with DD-PC1, PC-F1, PC-LA ,and NL3*. In general, PC-PK1 improves 
the overall agreement with data in comparison with the other interactions, especially for 
N = 82 and N = 126 isotones. The deviations are within 1 MeV for N = 82 isotones and 
2 MeV for both N = 20 and N = 50 isotones. A remarkable improvement in the binding 
energies as well as proper isospin dependence for N = 126 is found in the calculations with 
PC-PK1. In short, PC-PK1 provides better prediction for not only the binding energies but 
also its isospin dependence. 

B. Two-neutron separation energy 

From the binding energies, one can extract the two-neutron separation energies, S^n = 
Eb(N,Z) — Eb(N — 2,Z). In Fig. 4, the two-neutron separation energies for even-even 
O, Ca, Ni, and Sn isotopes predicted by PC-PK1 are shown in comparison with data [23] 
and those by DD-PC1, PC-F1, PC-LA, and NL3*. Generally speaking, similar as the 
other interactions, the calculation with PC-PK1 reproduces the experimental two-neutron 
separation energies quite well. For the oxygen isotopic chain, all the effective interactions 
predict the last bound neutron-rich nucleus as 28 contrary to experiment in which 24 is so 
far the last bound neutron-rich nucleus. One can also find that the deviations are large for 
the even-even 58_68 Ni which can be attributed to the underestimation of binding energies as 
shown in Fig. 2. Moreover, visible deviations between different predictions can be seen in 
the neutron-rich Sn isotopes, which requires future experimental confirmation. 

C. Single-particle level 

In Figs. 5, 6, the calculated single-particle energies for 16 0, 40 Ca, 132 Sn, and 208 Pb by 
PC-PK1 are shown in comparison with data [43] and those by DD-PC1, PC-F1, PC-LA, and 
NL3*. The experimental values are extracted from the single-nucleon separation energies or 
excitation energies [43]. The theoretical single-particle energies are the eigenvalues of the 
Dirac equation for nucleon. It should be kept in mind that the calculations are performed 
by neglecting particle- vibration coupling [44]. 

In Figs. 5, 6, it is clearly shown that the single-particle levels near the magic numbers and 
the corresponding shell gaps given by PC-PK1 are in good agreement with the experimental 
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values. In particular for 16 and 40 Ca, both the experimental proton and neutron single- 
particle spectra are well reproduced by PC-PK1. For 132 Sn and 208 Pb, the empirical levels 
close to the Fermi surface are also reproduced well. Moreover, the spurious shells at Z = 58 
( 132 Sn) and Z = 92 ( 208 Pb) are found for all the effective interactions, which may be improved 
by the inclusion of p-tensor couplings [45]. 

D. Charge radii and neutron skin thicknesses 

In Fig. 7, the charge radii for Sn and Pb isotopes predicted by PC-PK1 are shown in 
comparison with data [24, 25] and those by DD-PC1, PC-F1, PC-LA, and NL3*. It is seen 
that all the effective interactions reproduce the observed charge radii of Sn isotopes quite 
well (within 0.3%). For the Pb isotopes, the kink in the charge radii has been excellently 
reproduced by all the effective interactions. Quantitatively, the observed charge radii of Pb 
isotopes are reproduced by the calculations with DD-PC1, PC-F1, and NL3* within ~ 0.3%, 
while the calculations with PC-LA within ~ 0.5%. 

In Fig. 8, the neutron skin thicknesses for Sn isotopes and 208 Pb predicted by PC-PK1 
are shown in comparison with data [46, 47] and those by DD-PC1, PC-F1, PC-LA, and 
NL3*. For Sn isotopes, although PC-PK1 slightly overestimates the neutron skin thickness 
in comparison with DD-PC1, it nicely reproduces the isotopic trend. For 208 Pb, all the 
interactions except DD-PC1 give similar neutron skin thicknesses which are larger than the 
data deduced from antiprotonic atoms [48], polarized proton scattering [49, 50], elastic pro- 
ton scattering [51], proton-nucleus elastic scattering [47], and agree within the experimental 
error bar with that from inelastic a scattering [52] . The slightly overestimated neutron skin 
thicknesses are due to the enhanced symmetry energies for nuclear matter shown in Table IV. 
In overall, the DD-PC1 parametrization provides better description of experimental charge 
radii and neutron skin thickness due to its smaller symmetry energy at saturation density. 

VI. DEFORMED NUCLEI 

In this section, we will focus on the description of the binding energies and deformations 
for selected well-deformed even-even nuclei. In order to investigate the fission barrier, a 
constrained calculation is also carried out by taking 240 Pu as an example. 
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A. Binding energy and deformation 

The binding energies and quadrupole deformations of the ground states for Yb and U 
isotopes are investigated in axially deformed code with PC-PK1 in comparison with those 
with DD-PC1 and PC-F1. 

In the upper panels of Fig. 9, the deviations of the calculated binding energies with PC- 
PK1, DD-PC1, and PC-F1 from the data [23] are shown as circles, triangles, and squares 
respectively. 

Before taking into account the rotational correction for the binding energies, a system- 
atical underestimate of the binding energies around 3 MeV for both Yb and U isotopes 
is found for PC-PK1. For PC-F1, the difference between the calculated and the observed 
binding energy decreases monotonically with the isospin values, i.e., around 1 ~ —3 MeV 
for Yb isotopes and —2 ~ —5 MeV for U isotopes. As almost all the isotopes shown in Fig 9 
are used to adjust the parameters, the predicted binding energies by DD-PC1 are in good 
agreement with the data (within 1 MeV). 

After taking into account the energy correction due to the restoration of rotational sym- 
metry in the cranking approximation [53], the calculated results by PC-PK1 (filled circles) 
reproduces the data quite well for both Yb and U isotopes, and the deviations are within 1 
MeV. While the differences between the corrected binding energies given by PC-F1 (filled 
squares) and data are still large. Since DD-PC1 is adjusted to the binding energies of 64 
well-deformed nuclei, the rotational correction energy is not considered in the corresponding 
calculations. The energy correction due to the restoration of rotational symmetry can be 
taken into account with the microscopic treatment of angular momentum projection or the 
cranking approximation [54] . It is noted that difference between the cranking approximation 
and the angular momentum projection exists (for example, it reaches 1 MeV in 240 Pu) [55]. 
For simplicity of systematic calculations, only the cranking approximation is used here. 

In the lower panels of Fig. 9, the calculated quadrupole deformations for the ground 
states by PC-PK1, DD-PC1, and PC-F1 are given in comparison with the corresponding 
data [56]. It shows that the deformations and their corresponding evolutions with neutron 
number for both Yb and U isotopes are well reproduced by PC-PK1, DD-PC1, and PC-F1. 
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B. Fission barrier 



In Fig. 10, the potential energy curves for 240 Pu as functions of the quadrupole defor- 
mation p 2 are shown. The dashed and solid lines correspond to the axially-symmetric and 
the triaxial calculations with PC-PK1, respectively. In the case of triaxial calculation, the 
solid line refers to the minima for each (3 for the potential energy surface (PES) in the f3 — 7 
plane. For comparison, the axially-symmetric result given by PC-F1 is also included. 

It is found that the PC-PK1 provides not only a good description for the deformation 
of the ground state [56] but also the energy difference between the ground-state and the 
shape isomeric state [57]. Furthermore, after including the triaxiality, as shown in Fig 10, 
the fission barrier given by PC-PK1 is in agreement with the empirical value [58]. It should 
be noted that the pairing correlation plays an important role in the description of fission 
barrier. Discussion on the dependence of the fission barrier height on the pairing correlations 
can be found in Ref. [59]. 

VII. NUCLEAR EXCITED PROPERTIES 

As a test of the new parameter set PC-PK1 in the description of nuclear spectroscopic 
properties for low-lying excitation states, the collective excitation spectra and transition 
probabilities in 150 Nd as well as the characteristic collective observables for Nd isotopes will 
be calculated starting from a five-dimensional collective Hamiltonian in which the parameters 
are determined by constrained self-consistent RMF calculations for triaxial shapes [60-62]. 

In Fig. 11, the excitation energies and B(E2; L\ — > [L — 2}f) values for the yrast states 
in 150 Nd predicted by PC-PK1 are shown in comparison with data [63, 64] and those by DD- 
PC1 and PC-F1. It can be seen that all the effective interactions provide similar excitation 
energies and intraband B(E2) values for the yrast band and reproduce the data quite well. 

In Fig. 12, the characteristic collective observables R4/2 = E(4i)/E(2f) and B(E2; 2± — > 
0^) for Nd isotopes given by PC-PK1 are shown in comparison with data [63, 64] and those 
by DD-PC1 and PC-F1. It is found that all the parameter sets reproduce the data quite 
well. In particular, the calculations reproduce in detail the rapid increase of R4/2 and B(E2) 
with the neutron number, i.e., from R4/2 ~ 1-9 and B(E2) < 30 W.u. in near spherical 
144 Nd to R 4/2 ~ 3.3 and B(E2) > 150 W.u. in well-deformed 152 Nd. 
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It shows clearly that the new effective interaction PC-PK1 can provide a good description 
not only for the ground state properties in spherical and deformed nuclei but also for the 
nuclear spectroscopic properties of low-lying excitation states. 

VIII. SUMMARY 

In summary, a new parametrization PC-PK1 for the nuclear covariant energy density 
functional with nonlinear point-coupling interaction has been proposed by fitting to ob- 
servables of 60 selected spherical nuclei, including the binding energies, charge radii and 
empirical pairing gaps. By scaling the coupling constants in PC-PK1 in accordance with 
the QCD-based Lagrangian, it is found that all the nine parameters are natural. The success 
of PC-PK1 has been illustrated through the description for infinite nuclear matter and finite 
nuclei including the ground-state and low-lying excited states. 

For the spherical nuclei, PC-PK1 can provide better descriptions for the binding energies 
in comparison with DD-PC1, PC-F1, PC-LA, and NL3* sets. For neutron skin thicknesses, 
the DD-PC1 provides better description as compared with the other effective interactions 
due to its smaller symmetry energy at saturation density. 

Taking Yb and U isotopes as examples, it is found that the PC-PK1 reproduces the 
deformations and their corresponding evolutions with neutron number quite well. After 
taking into account the rotational correction energy in the cranking approximation, the 
binding energies given by PC-PK1 are in very good agreement with data within 1 MeV, 
which indicates that PC-PK1 achieves the same quality as DD-PC1 in the description for 
deformed nuclei. Moreover, PC-PK1 provides good description for isospin dependence of 
binding energy along either the isotopic or the isotonic chains, which makes it reliable for 
application in exotic nuclei. It is noted that the rotational correction energy evaluated using 
the cranking approximation may differ from that using angular momentum projection. 

Constrained calculations have also been performed for 240 Pu in order to investigate the 
fission barrier. It is found that the PC-PK1 provides not only a good description for the 
deformation of the ground state [56] but also the energy difference between the ground-state 
and the shape isomeric state [57]. Furthermore, after including the triaxiality, the fission 
barrier given by PC-PK1 is in agreement with the empirical value [58]. 

The predictive power of the PC-PK1 is also illustrated in the description for the collective 
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excitation spectra and transition probabilities in 150 Nd as well as the characteristic collec- 
tive observables for Nd isotopes in a five-dimensional collective Hamiltonian in which the 
parameters are determined by constrained calculations for triaxial shapes. There are also 
many extensions of nuclear covariant energy density functional theory beyond mean-field 
using projection techniques [65] and generator coordinate methods [66, 67]. More micro- 
scopic analysis of nuclear low-lying states in context of these frameworks with PC-PK1 is in 
progress. 

The density-dependent parametrization DD-PC1 is determined mainly from the masses 
of deformed nuclei and the EOS of nuclear matter. However, the calculations of the rear- 
rangement terms for the density-dependent parametrization can be nontrivial in some cases, 
in particular for RPA calculations. Here the nonlinear parametrization PC-PK1 has been 
optimized to the masses, charge radii and empirical pairing gaps for selected 60 spherical 
nuclei. It has been illustrated that the PC-PK1 can provide very good descriptions for both 
spherical and deformed nuclei. Therefore, the non-linear parametrization is very useful as 
it combines the simplicity with very good predictions for many nuclear properties. 
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TABLE I: The point-coupling constants and pairing strengths of PC-PK1 set. The corresponding 
QCD-scaled coupling constants q n are given in the last column as well. 



Coupling Constant 


Value 


Dimension 


Cln 


a s 


-3.96291 x 10" 4 


MeV" 2 


-1.695 


Ps 


8.6653 x 10- 11 


MeV~ 5 


1.628 


7s 


-3.80724 x KT 17 


MeV" 8 


-3.535 


Ss 


-1.09108 x 10~ 10 


MeV" 4 


-0.277 


ay 


2.6904 x 10" 4 


MeV" 2 


1.151 


IV 


-3.64219 x 10~ 18 


MeV" 8 


-0.338 


S v 


-4.32619 x 10~ 10 


MeV" 4 


-1.097 


aTv 


2.95018 x 10~ 5 


MeV" 2 


0.505 


Stv 


-4.11112 x 10~ 10 


MeV~ 4 


-4.171 


V n 


-349.5 


MeV fm 3 




v P 


-330 


MeV fm 3 





TABLE II: The calculated binding energies (in MeV) for selected spherical nuclei by PC- 
PK1 in comparison with the data [23] and those by DD-PC1 [15], PC-F1 [11], PC-LA 
[10], and NL3* [28]. The bold-faced quantities denote that the experimental values of 
the corresponding nuclei are used in the parametrization fitting. The root mean square 

N 



(rms) deviation A = y (E° xp — Ef 3 - 1 ) 2 /N and the root of relative square (rrs) deviation 



y — 1 % — are respectively listed in the last two rows. 



Nuclei 


Exp. 


PC-PK1 


DD-PC1 


PC-F1 


PC-LA 


NL3* 


16q 


127.619 


127.280 


128.527 


127.691 


127.407 


128.112 


18Q 


139.806 


140.223 


141.145 


140.028 


140.356 


140.504 


20Q 


151.370 


151.962 


152.790 


151.606 


152.228 


151.955 


22 Q 


162.026 


162.285 


163.141 


162.054 


162.665 


161.990 


18 Ne 


132.143 


132.088 


132.923 


132.216 


132.317 


132.494 



22 



TABLE II: (Continued). 



Nuclei 


Exp. 


PC-PK1 


DD-PC1 


PC-F1 


PC-LA 


NL3* 


20 Mg 


134.468 


134.563 


135.141 


134.613 


134.992 


134.786 


34 Si 


283.429 


284.727 


285.967 


285.067 


283.989 


283.236 


36 g 


308.714 


308.374 


309.305 


308.973 


307.221 


306.086 


38 Ar 


327.342 


327.107 


328.691 


328.540 


326.755 


325.379 


36 Ca 


281.360 


281.412 


281.878 


282.001 


280.454 


279.579 


38 Ca 


313.122 


313.230 


314.501 


314.415 


312.901 


311.669 


40 Ca 


342.052 


343.060 


345.113 


345.041 


343.202 


341.578 


42 Ca 


361.896 


363.142 


365.143 


364.411 


363.685 


361.547 


44 Ca 


380.960 


381.915 


383.967 


382.748 


382.789 


380.246 


46 Ca 


398.769 


399.451 


401.668 


400.060 


400.627 


397.718 


48 Ca 


415.990 


415.492 


417.973 


416.085 


416.969 


413.616 


50 Ca 


427.490 


426.937 


428.660 


427.302 


426.883 


424.445 


42 Ti 


346.905 


348.024 


349.848 


349.701 


348.626 


346.539 


50 Ti 


437.781 


436.445 


437.761 


436.171 


437.223 


434.389 


56 Ni 


483.992 


483.669 


481.447 


480.758 


481.826 


481.058 


58 Ni 


506.458 


503.636 


502.587 


501.646 


502.623 


501.342 


72 Ni 


613.169 


614.875 


617.071 


614.646 


614.486 


612.561 


84 Se 


727.343 


725.732 


728.792 


726.609 


727.605 


724.965 


86 Kr 


749.234 


747.939 


751.050 


749.427 


750.313 


747.055 


88 Sr 


768.468 


767.138 


770.240 


769.143 


769.742 


766.225 


90 Zr 


783.892 


783.033 


785.806 


785.348 


785.565 


782.336 


92 Mo 


796.508 


796.148 


798.308 


798.191 


798.719 


795.788 


94 Ru 


806.848 


807.034 


808.575 


808.731 


809.695 


807.019 


98 Cd 


821.067 


822.765 


823.162 


823.668 


825.580 


823.347 


100 Sn 


824.794 


827.715 


827.609 


828.156 


830.582 


828.529 


106 Sn 


893.868 


892.323 


893.469 


893.370 


895.447 


893.873 


108 Sn 


914.626 


913.179 


914.627 


914.236 


916.165 


914.665 
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TABLE II: (Continued). 



Nuclei 


Exp. 


PC-PK1 


DD-PC1 


PC-F1 


PC-LA 


NL3* 


112 Sn 


953.532 


951.831 


953.922 


953.367 


954.258 


952.866 


116 Sn 


988.684 


987.601 


990.019 


989.326 


989.016 


987.920 


120 Sn 


1020.546 


1020.415 


1022.902 


1021.704 


1020.767 


1020.014 


122 Sn 


1035.529 


1035.860 


1038.417 


1036.755 


1035.794 


1035.116 


124 Sn 


1049.963 


1050.715 


1053.402 


1051.160 


1050.327 


1049.631 


126 Sn 


1063.889 


1064.993 


1067.877 


1064.978 


1064.381 


1063.560 


128 Sn 


1077.346 


1078.688 


1081.835 


1078.234 


1077.945 


1076.885 


130 Sn 


1090.293 


1091.774 


1095.253 


1090.930 


1090.993 


1089.566 


132 Sn 


1102.851 


1104.202 


1108.096 


1103.057 


1103.484 


1101.551 


134 Sn 


1109.235 


1109.253 


1112.253 


1107.330 


1106.707 


1106.027 


i34 Te 


1123.434 


1124.205 


1128.176 


1124.193 


1124.613 


1122.859 


136 Xe 


1141.878 


1142.621 


1146.587 


1143.601 


1143.997 


1142.480 


138 Ba 


1158.292 


1159.381 


1163.283 


1161.245 


1161.575 


1160.331 


140 Ce 


1172.692 


1174.054 


1177.868 


1176.722 


1176.953 


1175.954 


142 Nd 


1185.141 


1185.938 


1189.537 


1189.138 


1189.292 


1188.002 


144 Sm 


1195.736 


1195.736 


1199.024 


1199.353 


1199.420 


1198.079 


146 Gd 


1204.435 


1203.712 


1206.614 


1207.635 


1207.687 


1206.449 


148 Dy 


1210.780 


1209.974 


1212.454 


1214.117 


1214.258 


1213.186 


i50 Er 


1215.331 


1214.624 


1216.686 


1218.943 


1219.236 


1218.343 


206 Hg 


1621.049 


1621.321 


1623.820 


1620.353 


1616.956 


1621.515 


200p b 


1576.354 


1574.885 


1577.817 


1575.666 


1575.769 


1578.189 


202p b 


1592.187 


1591.172 


1594.139 


1591.675 


1591.240 


1593.909 


204p b 


1607.506 


1607.068 


1610.026 


1607.325 


1606.187 


1609.199 


206p b 


1622.324 


1622.525 


1625.385 


1622.563 


1620.490 


1624.008 


208p b 


1636.430 


1637.438 


1640.008 


1637.241 


1633.865 


1638.237 


210p b 


1645.552 


1645.449 


1648.272 


1644.793 


1641.484 


1645.954 


212p b 


1654.514 


1653.425 


1656.428 


1652.275 


1648.887 


1653.546 
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TABLE II: (Continued). 



Nuclei 


Exp. 


PC-PK1 


DD-PC1 


PC-F1 


PC-LA 


NL3* 


214p b 


1663.291 


1661.397 


1664.481 


1659.697 


1656.073 


1661.056 


210p Q 




1 646 703 


1649.441 


1 647 760 

-L KJ^t I . i \J\J 


1 fi44 fi4S 


1 64S QQ5 


212 Rn 


1652 497 


1654 632 


1657 476 


1656 863 


1653 921 


1658 319 


214 Ra 


1658.315 


1661.172 


1664.092 


1664.512 


1661.709 


1666.174 


216 Xh 


1662.689 


1666.248 


1669.244 


1670.649 


1667.967 


1672.505 


21 8u 


1665.648 


1669.602 


1672.733 


1675.109 


1672.491 


1677.091 


A 




1.33 


3.09 


2.60 


2.64 


2.88 


5 




0.18% 


0.45% 


0.32% 


0.30% 


0.34% 




FIG. 1: (Color online) The binding energy per nucleon E/A for nuclear matter as a function of 
the baryon density p B given by PC-PK1, DD-PC1, PC-F1, PC-LA, NL3*, and PK1. The shaded 
area indicates the empirical value [34] and the filled diamonds present the microscopic results of 
the ab-initio variational calculation [38]. 
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TABLE III: The calculated charge radii (in fm) for selected spherical nuclei by PC-PK1 in com- 
parison with data [24, 25] and those by DD-PC1 [15], PC-F1 [11], PC-LA [10], and NL3* [28]. The 
bold-faced quantities denote that the experimental values of the corresponding nuclei are used in 

N 

the parametrization fitting. The root mean square (rms) deviation A = y (?^ xp — r? al ) 2 /iV and 

i 

N I ( r cxp - r ca1 ) 2 /(r cxp ) 2 
the root of relative square (rrs) deviation 6 = y — - l — - are respectively listed in 

i ' 

the last two rows. 



Nuclei 


Exp. 


PC-PK1 


DD-PC1 


PC-F1 


PC-LA 


NL3* 


16 


2.737 


2.7677 


2.7472 


2.7633 


2.7528 


2.7352 


40 Ca 


3.4852 


3.4815 


3.4566 


3.4777 


3.4678 


3.4704 


42 Ca 


3.5125 


3.4805 


3.4626 


3.4778 


3.4729 


3.4672 


44 Ca 


3.5231 


3.4826 


3.4709 


3.4809 


3.4810 


3.4672 


46 Ca 


3.5022 


3.4865 


3.4806 


3.4860 


3.4912 


3.4693 


48 Ca 


3.4837 


3.4890 


3.4895 


3.4906 


3.5023 


3.4705 


50 Ti 


3.573 


3.5558 


3.5696 


3.5664 


3.5868 


3.5442 


58 M 


3.7827 


3.7372 


3.7761 


3.7645 


3.8065 


3.7399 


88 Sr 


4.2036 


4.2247 


4.2231 


4.2269 


4.2379 


4.2159 


90 Zr 


4.2720 


4.2695 


4.2664 


4.2724 


4.2847 


4.2636 


92 Mo 


4.3170 


4.3125 


4.3140 


4.3192 


4.3333 


4.3087 


112 Sn 


4.5957 


4.5801 


4.5894 


4.5870 


4.6044 


4.5753 


116 Sn 


4.6257 


4.6121 


4.6174 


4.6168 


4.6307 


4.6039 


122 Sn 


4.6633 


4.6561 


4.6579 


4.6549 


4.6728 


4.6430 


124 Sn 


4.6739 


4.6694 


4.6714 


4.6677 


4.6864 


4.6554 


138 Ba 


4.8348 


4.8508 


4.8511 


4.8494 


4.8667 


4.8369 


140 Ce 


4.8774 


4.8879 


4.8879 


4.8871 


4.9037 


4.8748 


144 Sm 


4.9525 


4.9544 


4.9521 


4.9547 


4.9676 


4.9484 


202p b 


5.4772 


5.4908 


5.4869 


5.4892 


5.4996 


5.4825 


204p b 


5.4861 


5.5005 


5.4962 


5.4987 


5.5112 


5.4916 


206p b 


5.4946 


5.5098 


5.5049 


5.5078 


5.5200 


5.5004 


208p b 


5.5046 


5.5185 


5.5129 


5.5162 


5.5279 


5.5087 


214p b 


5.5622 


5.5798 


5.5711 


5.5762 


5.5813 


5.5699 


A 




0.019 


G>§19 


0.017 


0.023 


0.022 


6 




0.53% 


0.51% 


0.45% 


0.55% 


0.60% 




FIG. 2: (Color online) Deviations of the calculated binding energies for Ca, Ni, Sn, and Pb isotopes 
by PC-PK1 from the data [23] in comparison with those by DD-PC1, PC-F1, PC-LA, and NL3*. 




FIG. 3: (Color online) Same as Fig. 2 but for the N=20, N=50, N=82 and N=126 isotones. 
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TABLE IV: The predicted saturation properties for nuclear matter by PC-PK1 in comparison with 
those by DD-PC1, PC-F1, PC-LA, NL3*, and PK1. 





PP PK1 


r>n ppi 


pr 1 fi 

I V. — I 1 


PC T A 


NT "3* 

IN IjO 


PK1 


rirx ( "fm ^ ^ 
PO \ IU1 ) 




U. Ui 


n i p.i 

U. -LUX 


n 1 48 

U. ±*-kO 


U. 1JU 


n i4s 

U. L<-tO 


S2j I -til ivit; v f 


1 1 9 


1 fi Ofi 


1 fi 1 7 


1 fi 1 "3 


1 fi ^1 


1 fi 97 


M*/M 


0.59 


0.58 


0.61 


0.58 


0.59 


0.60 


M* L /M 


0.65 


0.64 


0.67 


0.64 


0.65 


0.66 


K (MeV) 


238 


230 


255 


264 


258 


283 


£ sym (MeV) 


35.6 


33 


37.8 


37.2 


38.7 


37.6 


L(MeV) 


113 


70 


117 


108 


123 


116 


^asy(MeV) 


-583 


-528 


-627 


-709 


-630 


-641 
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FIG. 4: (Color online) The calculated two- neutron separation energies for O, Ca, Ni, and Sn 
isotopes by PC-PK1 in comparison with data [23] and those by DD-PC1, PC-F1, PC-LA, and 
NL3*. 
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FIG. 5: (Color online) The calculated single-particle energies for 16 and 40 Ca by PC-PK1 in 
comparison with data [43] and those by DD-PC1, PC-F1, PC-LA, and NL3*. 
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FIG. 6: (Color online) Same as Fig. 5 but for 132 Sn and 208 Pb. 
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FIG. 7: (Color online) The calculated charge radii for Sn and Pb isotopes by the PC-PK1 in 
comparison with data [24, 25] and those by DD-PC1, PC-F1, PC-LA, and NL3*. 
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FIG. 8: (Color online) The calculated neutron skin thicknesses for Sn isotopes and 208 Pb by 
PC-PK1 in comparison with data [46, 47] and those by DD-PC1, PC-F1, PC-LA, and NL3*. 
In the lower panel, the data for 208 Pb deduced from antiprotonic atoms [48], polarized proton 
scattering [49, 50], elastic proton scattering [51], inelastic a scattering [52], and proton- nucleus 
elastic scattering [47] are shown from left to right respectively. 
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FIG. 9: (Color online) Deviations of the calculated binding energies from the data [23] for Yb and 
U isotopes in axially deformed code by PC-PK1, DD-PC1 and PC-F1 (upper panel) as well as the 
corresponding calculated ground-state deformations in comparison with data [56] (lower panel). 
The filled circles and squares in the upper panels correspond to the rotational corrected ones given 
by PC-PK1 and PC-F1 respectively. 
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FIG. 10: (Color online) The potential energy curves for 240 Pu as functions of the quadrupole 
deformation fa in the calculations with PC-PK1. The dashed and solid lines correspond to the 
axial results and the triaxial results, i.e., the minima for each j3 for the potential energy surface 
(PES) in the /3 — 7 plane, respectively. For comparison, the axially-symmetric result given by 
PC-F1 is also included as dot-dashed line. The data for the ground-state deformation [56], the 
barrier height [58], and the energy of the fission isomer [57] are respectively indicated by an arrow, 
a diamond, and a square. To guide the eyes, the diamond and square are respectively set at 
fa = 0.64 and fa = 0.95. 
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FIG. 11: (Color online) The predicted excitation energies (upper panel) and B(E2; Lf — > [L — 2]f) 
values (lower panel) for the yrast states in 150 Nd by PC-PK1 in comparison with data [63, 64] 
and those by DD-PC1 and PC-F1. The energies and B(E2) values are respectively normalized to 
E(2+) and B(E2;2+ -> 0+). 
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FIG. 12: (Color online) The predicted characteristic collective observables R4/2 = E(AJ)/E(2J) 
and B(E2;2f — >• Of) (in W.u.) for Nd isotopes by PC-PK1 in comparison with data [63, 64] and 
those by DD-PC1 and PC-F1. 
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